3.116 \(\int \frac{x^3}{\log ^3(c (a+b x^2)^p)} \, dx\)

Optimal. Leaf size=204 \[ \frac{\left (a+b x^2\right )^2 \left (c \left (a+b x^2\right )^p\right )^{-2/p} \text{Ei}\left (\frac{2 \log \left (c \left (b x^2+a\right )^p\right )}{p}\right )}{b^2 p^3}-\frac{a \left (a+b x^2\right ) \left (c \left (a+b x^2\right )^p\right )^{-1/p} \text{Ei}\left (\frac{\log \left (c \left (b x^2+a\right )^p\right )}{p}\right )}{4 b^2 p^3}-\frac{a \left (a+b x^2\right )}{4 b^2 p^2 \log \left (c \left (a+b x^2\right )^p\right )}-\frac{x^2 \left (a+b x^2\right )}{2 b p^2 \log \left (c \left (a+b x^2\right )^p\right )}-\frac{x^2 \left (a+b x^2\right )}{4 b p \log ^2\left (c \left (a+b x^2\right )^p\right )} \]

[Out]

-(a*(a + b*x^2)*ExpIntegralEi[Log[c*(a + b*x^2)^p]/p])/(4*b^2*p^3*(c*(a + b*x^2)^p)^p^(-1)) + ((a + b*x^2)^2*E
xpIntegralEi[(2*Log[c*(a + b*x^2)^p])/p])/(b^2*p^3*(c*(a + b*x^2)^p)^(2/p)) - (x^2*(a + b*x^2))/(4*b*p*Log[c*(
a + b*x^2)^p]^2) - (a*(a + b*x^2))/(4*b^2*p^2*Log[c*(a + b*x^2)^p]) - (x^2*(a + b*x^2))/(2*b*p^2*Log[c*(a + b*
x^2)^p])

________________________________________________________________________________________

Rubi [A]  time = 0.288515, antiderivative size = 204, normalized size of antiderivative = 1., number of steps used = 18, number of rules used = 9, integrand size = 18, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.5, Rules used = {2454, 2400, 2399, 2389, 2300, 2178, 2390, 2310, 2297} \[ \frac{\left (a+b x^2\right )^2 \left (c \left (a+b x^2\right )^p\right )^{-2/p} \text{Ei}\left (\frac{2 \log \left (c \left (b x^2+a\right )^p\right )}{p}\right )}{b^2 p^3}-\frac{a \left (a+b x^2\right ) \left (c \left (a+b x^2\right )^p\right )^{-1/p} \text{Ei}\left (\frac{\log \left (c \left (b x^2+a\right )^p\right )}{p}\right )}{4 b^2 p^3}-\frac{a \left (a+b x^2\right )}{4 b^2 p^2 \log \left (c \left (a+b x^2\right )^p\right )}-\frac{x^2 \left (a+b x^2\right )}{2 b p^2 \log \left (c \left (a+b x^2\right )^p\right )}-\frac{x^2 \left (a+b x^2\right )}{4 b p \log ^2\left (c \left (a+b x^2\right )^p\right )} \]

Antiderivative was successfully verified.

[In]

Int[x^3/Log[c*(a + b*x^2)^p]^3,x]

[Out]

-(a*(a + b*x^2)*ExpIntegralEi[Log[c*(a + b*x^2)^p]/p])/(4*b^2*p^3*(c*(a + b*x^2)^p)^p^(-1)) + ((a + b*x^2)^2*E
xpIntegralEi[(2*Log[c*(a + b*x^2)^p])/p])/(b^2*p^3*(c*(a + b*x^2)^p)^(2/p)) - (x^2*(a + b*x^2))/(4*b*p*Log[c*(
a + b*x^2)^p]^2) - (a*(a + b*x^2))/(4*b^2*p^2*Log[c*(a + b*x^2)^p]) - (x^2*(a + b*x^2))/(2*b*p^2*Log[c*(a + b*
x^2)^p])

Rule 2454

Int[((a_.) + Log[(c_.)*((d_) + (e_.)*(x_)^(n_))^(p_.)]*(b_.))^(q_.)*(x_)^(m_.), x_Symbol] :> Dist[1/n, Subst[I
nt[x^(Simplify[(m + 1)/n] - 1)*(a + b*Log[c*(d + e*x)^p])^q, x], x, x^n], x] /; FreeQ[{a, b, c, d, e, m, n, p,
q}, x] && IntegerQ[Simplify[(m + 1)/n]] && (GtQ[(m + 1)/n, 0] || IGtQ[q, 0]) &&  !(EqQ[q, 1] && ILtQ[n, 0] &&
IGtQ[m, 0])

Rule 2400

Int[((a_.) + Log[(c_.)*((d_) + (e_.)*(x_))^(n_.)]*(b_.))^(p_)*((f_.) + (g_.)*(x_))^(q_.), x_Symbol] :> Simp[((
d + e*x)*(f + g*x)^q*(a + b*Log[c*(d + e*x)^n])^(p + 1))/(b*e*n*(p + 1)), x] + (-Dist[(q + 1)/(b*n*(p + 1)), I
nt[(f + g*x)^q*(a + b*Log[c*(d + e*x)^n])^(p + 1), x], x] + Dist[(q*(e*f - d*g))/(b*e*n*(p + 1)), Int[(f + g*x
)^(q - 1)*(a + b*Log[c*(d + e*x)^n])^(p + 1), x], x]) /; FreeQ[{a, b, c, d, e, f, g, n}, x] && NeQ[e*f - d*g,
0] && LtQ[p, -1] && GtQ[q, 0]

Rule 2399

Int[((f_.) + (g_.)*(x_))^(q_.)/((a_.) + Log[(c_.)*((d_) + (e_.)*(x_))^(n_.)]*(b_.)), x_Symbol] :> Int[ExpandIn
tegrand[(f + g*x)^q/(a + b*Log[c*(d + e*x)^n]), x], x] /; FreeQ[{a, b, c, d, e, f, g, n}, x] && NeQ[e*f - d*g,
0] && IGtQ[q, 0]

Rule 2389

Int[((a_.) + Log[(c_.)*((d_) + (e_.)*(x_))^(n_.)]*(b_.))^(p_.), x_Symbol] :> Dist[1/e, Subst[Int[(a + b*Log[c*
x^n])^p, x], x, d + e*x], x] /; FreeQ[{a, b, c, d, e, n, p}, x]

Rule 2300

Int[((a_.) + Log[(c_.)*(x_)^(n_.)]*(b_.))^(p_), x_Symbol] :> Dist[x/(n*(c*x^n)^(1/n)), Subst[Int[E^(x/n)*(a +
b*x)^p, x], x, Log[c*x^n]], x] /; FreeQ[{a, b, c, n, p}, x]

Rule 2178

Int[(F_)^((g_.)*((e_.) + (f_.)*(x_)))/((c_.) + (d_.)*(x_)), x_Symbol] :> Simp[(F^(g*(e - (c*f)/d))*ExpIntegral
Ei[(f*g*(c + d*x)*Log[F])/d])/d, x] /; FreeQ[{F, c, d, e, f, g}, x] &&  !$UseGamma === True

Rule 2390

Int[((a_.) + Log[(c_.)*((d_) + (e_.)*(x_))^(n_.)]*(b_.))^(p_.)*((f_) + (g_.)*(x_))^(q_.), x_Symbol] :> Dist[1/
e, Subst[Int[((f*x)/d)^q*(a + b*Log[c*x^n])^p, x], x, d + e*x], x] /; FreeQ[{a, b, c, d, e, f, g, n, p, q}, x]
&& EqQ[e*f - d*g, 0]

Rule 2310

Int[((a_.) + Log[(c_.)*(x_)^(n_.)]*(b_.))^(p_)*((d_.)*(x_))^(m_.), x_Symbol] :> Dist[(d*x)^(m + 1)/(d*n*(c*x^n
)^((m + 1)/n)), Subst[Int[E^(((m + 1)*x)/n)*(a + b*x)^p, x], x, Log[c*x^n]], x] /; FreeQ[{a, b, c, d, m, n, p}
, x]

Rule 2297

Int[((a_.) + Log[(c_.)*(x_)^(n_.)]*(b_.))^(p_), x_Symbol] :> Simp[(x*(a + b*Log[c*x^n])^(p + 1))/(b*n*(p + 1))
, x] - Dist[1/(b*n*(p + 1)), Int[(a + b*Log[c*x^n])^(p + 1), x], x] /; FreeQ[{a, b, c, n}, x] && LtQ[p, -1] &&
IntegerQ[2*p]

Rubi steps

\begin{align*} \int \frac{x^3}{\log ^3\left (c \left (a+b x^2\right )^p\right )} \, dx &=\frac{1}{2} \operatorname{Subst}\left (\int \frac{x}{\log ^3\left (c (a+b x)^p\right )} \, dx,x,x^2\right )\\ &=-\frac{x^2 \left (a+b x^2\right )}{4 b p \log ^2\left (c \left (a+b x^2\right )^p\right )}+\frac{\operatorname{Subst}\left (\int \frac{x}{\log ^2\left (c (a+b x)^p\right )} \, dx,x,x^2\right )}{2 p}+\frac{a \operatorname{Subst}\left (\int \frac{1}{\log ^2\left (c (a+b x)^p\right )} \, dx,x,x^2\right )}{4 b p}\\ &=-\frac{x^2 \left (a+b x^2\right )}{4 b p \log ^2\left (c \left (a+b x^2\right )^p\right )}-\frac{x^2 \left (a+b x^2\right )}{2 b p^2 \log \left (c \left (a+b x^2\right )^p\right )}+\frac{\operatorname{Subst}\left (\int \frac{x}{\log \left (c (a+b x)^p\right )} \, dx,x,x^2\right )}{p^2}+\frac{a \operatorname{Subst}\left (\int \frac{1}{\log \left (c (a+b x)^p\right )} \, dx,x,x^2\right )}{2 b p^2}+\frac{a \operatorname{Subst}\left (\int \frac{1}{\log ^2\left (c x^p\right )} \, dx,x,a+b x^2\right )}{4 b^2 p}\\ &=-\frac{x^2 \left (a+b x^2\right )}{4 b p \log ^2\left (c \left (a+b x^2\right )^p\right )}-\frac{a \left (a+b x^2\right )}{4 b^2 p^2 \log \left (c \left (a+b x^2\right )^p\right )}-\frac{x^2 \left (a+b x^2\right )}{2 b p^2 \log \left (c \left (a+b x^2\right )^p\right )}+\frac{\operatorname{Subst}\left (\int \left (-\frac{a}{b \log \left (c (a+b x)^p\right )}+\frac{a+b x}{b \log \left (c (a+b x)^p\right )}\right ) \, dx,x,x^2\right )}{p^2}+\frac{a \operatorname{Subst}\left (\int \frac{1}{\log \left (c x^p\right )} \, dx,x,a+b x^2\right )}{4 b^2 p^2}+\frac{a \operatorname{Subst}\left (\int \frac{1}{\log \left (c x^p\right )} \, dx,x,a+b x^2\right )}{2 b^2 p^2}\\ &=-\frac{x^2 \left (a+b x^2\right )}{4 b p \log ^2\left (c \left (a+b x^2\right )^p\right )}-\frac{a \left (a+b x^2\right )}{4 b^2 p^2 \log \left (c \left (a+b x^2\right )^p\right )}-\frac{x^2 \left (a+b x^2\right )}{2 b p^2 \log \left (c \left (a+b x^2\right )^p\right )}+\frac{\operatorname{Subst}\left (\int \frac{a+b x}{\log \left (c (a+b x)^p\right )} \, dx,x,x^2\right )}{b p^2}-\frac{a \operatorname{Subst}\left (\int \frac{1}{\log \left (c (a+b x)^p\right )} \, dx,x,x^2\right )}{b p^2}+\frac{\left (a \left (a+b x^2\right ) \left (c \left (a+b x^2\right )^p\right )^{-1/p}\right ) \operatorname{Subst}\left (\int \frac{e^{\frac{x}{p}}}{x} \, dx,x,\log \left (c \left (a+b x^2\right )^p\right )\right )}{4 b^2 p^3}+\frac{\left (a \left (a+b x^2\right ) \left (c \left (a+b x^2\right )^p\right )^{-1/p}\right ) \operatorname{Subst}\left (\int \frac{e^{\frac{x}{p}}}{x} \, dx,x,\log \left (c \left (a+b x^2\right )^p\right )\right )}{2 b^2 p^3}\\ &=\frac{3 a \left (a+b x^2\right ) \left (c \left (a+b x^2\right )^p\right )^{-1/p} \text{Ei}\left (\frac{\log \left (c \left (a+b x^2\right )^p\right )}{p}\right )}{4 b^2 p^3}-\frac{x^2 \left (a+b x^2\right )}{4 b p \log ^2\left (c \left (a+b x^2\right )^p\right )}-\frac{a \left (a+b x^2\right )}{4 b^2 p^2 \log \left (c \left (a+b x^2\right )^p\right )}-\frac{x^2 \left (a+b x^2\right )}{2 b p^2 \log \left (c \left (a+b x^2\right )^p\right )}+\frac{\operatorname{Subst}\left (\int \frac{x}{\log \left (c x^p\right )} \, dx,x,a+b x^2\right )}{b^2 p^2}-\frac{a \operatorname{Subst}\left (\int \frac{1}{\log \left (c x^p\right )} \, dx,x,a+b x^2\right )}{b^2 p^2}\\ &=\frac{3 a \left (a+b x^2\right ) \left (c \left (a+b x^2\right )^p\right )^{-1/p} \text{Ei}\left (\frac{\log \left (c \left (a+b x^2\right )^p\right )}{p}\right )}{4 b^2 p^3}-\frac{x^2 \left (a+b x^2\right )}{4 b p \log ^2\left (c \left (a+b x^2\right )^p\right )}-\frac{a \left (a+b x^2\right )}{4 b^2 p^2 \log \left (c \left (a+b x^2\right )^p\right )}-\frac{x^2 \left (a+b x^2\right )}{2 b p^2 \log \left (c \left (a+b x^2\right )^p\right )}+\frac{\left (\left (a+b x^2\right )^2 \left (c \left (a+b x^2\right )^p\right )^{-2/p}\right ) \operatorname{Subst}\left (\int \frac{e^{\frac{2 x}{p}}}{x} \, dx,x,\log \left (c \left (a+b x^2\right )^p\right )\right )}{b^2 p^3}-\frac{\left (a \left (a+b x^2\right ) \left (c \left (a+b x^2\right )^p\right )^{-1/p}\right ) \operatorname{Subst}\left (\int \frac{e^{\frac{x}{p}}}{x} \, dx,x,\log \left (c \left (a+b x^2\right )^p\right )\right )}{b^2 p^3}\\ &=-\frac{a \left (a+b x^2\right ) \left (c \left (a+b x^2\right )^p\right )^{-1/p} \text{Ei}\left (\frac{\log \left (c \left (a+b x^2\right )^p\right )}{p}\right )}{4 b^2 p^3}+\frac{\left (a+b x^2\right )^2 \left (c \left (a+b x^2\right )^p\right )^{-2/p} \text{Ei}\left (\frac{2 \log \left (c \left (a+b x^2\right )^p\right )}{p}\right )}{b^2 p^3}-\frac{x^2 \left (a+b x^2\right )}{4 b p \log ^2\left (c \left (a+b x^2\right )^p\right )}-\frac{a \left (a+b x^2\right )}{4 b^2 p^2 \log \left (c \left (a+b x^2\right )^p\right )}-\frac{x^2 \left (a+b x^2\right )}{2 b p^2 \log \left (c \left (a+b x^2\right )^p\right )}\\ \end{align*}

Mathematica [A]  time = 0.193008, size = 185, normalized size = 0.91 \[ -\frac{\left (a+b x^2\right ) \left (c \left (a+b x^2\right )^p\right )^{-2/p} \left (a \left (c \left (a+b x^2\right )^p\right )^{\frac{1}{p}} \log ^2\left (c \left (a+b x^2\right )^p\right ) \text{Ei}\left (\frac{\log \left (c \left (b x^2+a\right )^p\right )}{p}\right )-4 \left (a+b x^2\right ) \log ^2\left (c \left (a+b x^2\right )^p\right ) \text{Ei}\left (\frac{2 \log \left (c \left (b x^2+a\right )^p\right )}{p}\right )+p \left (c \left (a+b x^2\right )^p\right )^{2/p} \left (\left (a+2 b x^2\right ) \log \left (c \left (a+b x^2\right )^p\right )+b p x^2\right )\right )}{4 b^2 p^3 \log ^2\left (c \left (a+b x^2\right )^p\right )} \]

Antiderivative was successfully verified.

[In]

Integrate[x^3/Log[c*(a + b*x^2)^p]^3,x]

[Out]

-((a + b*x^2)*(a*(c*(a + b*x^2)^p)^p^(-1)*ExpIntegralEi[Log[c*(a + b*x^2)^p]/p]*Log[c*(a + b*x^2)^p]^2 - 4*(a
+ b*x^2)*ExpIntegralEi[(2*Log[c*(a + b*x^2)^p])/p]*Log[c*(a + b*x^2)^p]^2 + p*(c*(a + b*x^2)^p)^(2/p)*(b*p*x^2
+ (a + 2*b*x^2)*Log[c*(a + b*x^2)^p])))/(4*b^2*p^3*(c*(a + b*x^2)^p)^(2/p)*Log[c*(a + b*x^2)^p]^2)

________________________________________________________________________________________

Maple [F]  time = 5.264, size = 0, normalized size = 0. \begin{align*} \int{\frac{{x}^{3}}{ \left ( \ln \left ( c \left ( b{x}^{2}+a \right ) ^{p} \right ) \right ) ^{3}}}\, dx \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x^3/ln(c*(b*x^2+a)^p)^3,x)

[Out]

int(x^3/ln(c*(b*x^2+a)^p)^3,x)

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} -\frac{b^{2}{\left (p + 2 \, \log \left (c\right )\right )} x^{4} + a b{\left (p + 3 \, \log \left (c\right )\right )} x^{2} + a^{2} \log \left (c\right ) +{\left (2 \, b^{2} x^{4} + 3 \, a b x^{2} + a^{2}\right )} \log \left ({\left (b x^{2} + a\right )}^{p}\right )}{4 \,{\left (b^{2} p^{2} \log \left ({\left (b x^{2} + a\right )}^{p}\right )^{2} + 2 \, b^{2} p^{2} \log \left ({\left (b x^{2} + a\right )}^{p}\right ) \log \left (c\right ) + b^{2} p^{2} \log \left (c\right )^{2}\right )}} + \int \frac{4 \, b x^{3} + 3 \, a x}{2 \,{\left (b p^{2} \log \left ({\left (b x^{2} + a\right )}^{p}\right ) + b p^{2} \log \left (c\right )\right )}}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^3/log(c*(b*x^2+a)^p)^3,x, algorithm="maxima")

[Out]

-1/4*(b^2*(p + 2*log(c))*x^4 + a*b*(p + 3*log(c))*x^2 + a^2*log(c) + (2*b^2*x^4 + 3*a*b*x^2 + a^2)*log((b*x^2
+ a)^p))/(b^2*p^2*log((b*x^2 + a)^p)^2 + 2*b^2*p^2*log((b*x^2 + a)^p)*log(c) + b^2*p^2*log(c)^2) + integrate(1
/2*(4*b*x^3 + 3*a*x)/(b*p^2*log((b*x^2 + a)^p) + b*p^2*log(c)), x)

________________________________________________________________________________________

Fricas [A]  time = 2.20534, size = 629, normalized size = 3.08 \begin{align*} -\frac{{\left (a p^{2} \log \left (b x^{2} + a\right )^{2} + 2 \, a p \log \left (b x^{2} + a\right ) \log \left (c\right ) + a \log \left (c\right )^{2}\right )} c^{\left (\frac{1}{p}\right )} \logintegral \left ({\left (b x^{2} + a\right )} c^{\left (\frac{1}{p}\right )}\right ) +{\left (b^{2} p^{2} x^{4} + a b p^{2} x^{2} +{\left (2 \, b^{2} p^{2} x^{4} + 3 \, a b p^{2} x^{2} + a^{2} p^{2}\right )} \log \left (b x^{2} + a\right ) +{\left (2 \, b^{2} p x^{4} + 3 \, a b p x^{2} + a^{2} p\right )} \log \left (c\right )\right )} c^{\frac{2}{p}} - 4 \,{\left (p^{2} \log \left (b x^{2} + a\right )^{2} + 2 \, p \log \left (b x^{2} + a\right ) \log \left (c\right ) + \log \left (c\right )^{2}\right )} \logintegral \left ({\left (b^{2} x^{4} + 2 \, a b x^{2} + a^{2}\right )} c^{\frac{2}{p}}\right )}{4 \,{\left (b^{2} p^{5} \log \left (b x^{2} + a\right )^{2} + 2 \, b^{2} p^{4} \log \left (b x^{2} + a\right ) \log \left (c\right ) + b^{2} p^{3} \log \left (c\right )^{2}\right )} c^{\frac{2}{p}}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^3/log(c*(b*x^2+a)^p)^3,x, algorithm="fricas")

[Out]

-1/4*((a*p^2*log(b*x^2 + a)^2 + 2*a*p*log(b*x^2 + a)*log(c) + a*log(c)^2)*c^(1/p)*log_integral((b*x^2 + a)*c^(
1/p)) + (b^2*p^2*x^4 + a*b*p^2*x^2 + (2*b^2*p^2*x^4 + 3*a*b*p^2*x^2 + a^2*p^2)*log(b*x^2 + a) + (2*b^2*p*x^4 +
3*a*b*p*x^2 + a^2*p)*log(c))*c^(2/p) - 4*(p^2*log(b*x^2 + a)^2 + 2*p*log(b*x^2 + a)*log(c) + log(c)^2)*log_in
tegral((b^2*x^4 + 2*a*b*x^2 + a^2)*c^(2/p)))/((b^2*p^5*log(b*x^2 + a)^2 + 2*b^2*p^4*log(b*x^2 + a)*log(c) + b^
2*p^3*log(c)^2)*c^(2/p))

________________________________________________________________________________________

Sympy [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int \frac{x^{3}}{\log{\left (c \left (a + b x^{2}\right )^{p} \right )}^{3}}\, dx \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x**3/ln(c*(b*x**2+a)**p)**3,x)

[Out]

Integral(x**3/log(c*(a + b*x**2)**p)**3, x)

________________________________________________________________________________________

Giac [B]  time = 1.24803, size = 1134, normalized size = 5.56 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^3/log(c*(b*x^2+a)^p)^3,x, algorithm="giac")

[Out]

-1/4*(2*(b*x^2 + a)^2*p^2*log(b*x^2 + a)/(b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*log(c) + b*p^3*log(c
)^2) - (b*x^2 + a)*a*p^2*log(b*x^2 + a)/(b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*log(c) + b*p^3*log(c)
^2) + a*p^2*Ei(log(c)/p + log(b*x^2 + a))*log(b*x^2 + a)^2/((b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*l
og(c) + b*p^3*log(c)^2)*c^(1/p)) + (b*x^2 + a)^2*p^2/(b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*log(c) +
b*p^3*log(c)^2) - (b*x^2 + a)*a*p^2/(b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*log(c) + b*p^3*log(c)^2)
- 4*p^2*Ei(2*log(c)/p + 2*log(b*x^2 + a))*log(b*x^2 + a)^2/((b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*
log(c) + b*p^3*log(c)^2)*c^(2/p)) + 2*(b*x^2 + a)^2*p*log(c)/(b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*
log(c) + b*p^3*log(c)^2) - (b*x^2 + a)*a*p*log(c)/(b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*log(c) + b*
p^3*log(c)^2) + 2*a*p*Ei(log(c)/p + log(b*x^2 + a))*log(b*x^2 + a)*log(c)/((b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*l
og(b*x^2 + a)*log(c) + b*p^3*log(c)^2)*c^(1/p)) - 8*p*Ei(2*log(c)/p + 2*log(b*x^2 + a))*log(b*x^2 + a)*log(c)/
((b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*log(c) + b*p^3*log(c)^2)*c^(2/p)) + a*Ei(log(c)/p + log(b*x^
2 + a))*log(c)^2/((b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*log(c) + b*p^3*log(c)^2)*c^(1/p)) - 4*Ei(2*
log(c)/p + 2*log(b*x^2 + a))*log(c)^2/((b*p^5*log(b*x^2 + a)^2 + 2*b*p^4*log(b*x^2 + a)*log(c) + b*p^3*log(c)^
2)*c^(2/p)))/b